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Chapter  V  -  A  Gas 'illne  Blond  Problem^ 

Alan  S.  Manne 

1.  Basic  econoadc  aasumptiona 

The  precelln^-'  chapter  has  dealt  »rith  one  type  of  Intra-refinery  ^conoxlc 
balance  problem  —  the  operation  of  a  typical  kind  of  conversion 
present  one  concentrates  upon  another  such  intra- re  finery  study  —  the 
selection  of  an  optimal  (gasoline  Meniinp  scheduleyj  Kot  until  tJy  following 
chapter  will  any  attempt  oe  mjido  to  liscuss  the  operations  of  an  inte>Qatei 
refinery. 

There  were  two  primary  reasons  for  sinp.lin,;  out  the  blond In^^  problems  of 
the  Union  Oil  Company  of  California:  First,  this  operatior,  is  itself  a  complex 
larg^e-scale  one  --  p:rossln^  an  amount  of  the  omer  of  $100  mllll)r.s  annually. 
And  second,  the  data  available  for  this  analysis  were  of  tolerable  precision.. 

It  was  a>^reed  that  there  is  little  point  In  rraKinp;  an  elaborate  econondc  study 
when  the  basic  numbers  are  shaky. 

TOf'ether  with  wmbers  of  the  Manufacturin':  Economics  Division  of  Union  Oil, 
the  problem  was  formulated  in  the  followin’  terms  (see  the  schematic  lla^ram, 
Fi^re  10):  First,  t^e  quantities  and  qualities  of  1^  raw  .^solino  stocks 
Ware  aasuised  given.  Together  with  tetraethyl  l“ad  (a  piirchas'*  1  item),  these 
serve  as  the  inputs  to  the  gasoline  bleniing  operation.  Tlie  ^nput  materials 
are  coablnei  to  form  three  basic  types  of  products:  7bX,  premium  grade  gaso¬ 
line;  76,  regular  grade  gasoline;  and  fuel  oil  cutter  stock.  As  a  first 
approximation,  it  was  assuraed  that  each  of  these  products  coul  i  be  solo  at 
refinery  roalliat ions  that  were  independent  of  the  total  quantity  T*rketed. 

The  scheduling  problem  consists  of  assigning  each  raw  gasolire  stock  among 

1.  In  connection  with  tUs  section,  a  particular  lebt  of  gratitude  is  o’<#ed  to 
Messrs.  Reaugh,  McCreery,  and  Norton  of  the  Manufacturin.’,  Economics  Division 
of  the  Union  Oil  Company  of  California,  Without  their  encouragement  and  a-tlve 
help,  it  would  net  have  oeen  possible  to  carry  through  the  work  described  here. 


the  three  potentlaJ.  products,  aud  al8<j  setting  tfu?  tetraethyl  lesd  levels 
for  the  tvo  actor  fuels.  This  assignment  has  to  take  Into  scccwnt  not  only 
the  ebcolutc  ^lantlties  (.f  materials  that  arc  available,  but  sleo  certain 
product  cpeciflcatioiis  on  Liie  two  finished  gstoiineB.’ 

lach  of  t.ucse  aojumptlons  musL  be  reviewed  with  care.  Consider,  for 
oxa-aplo,  the  antter  of  blending  stocks.  With  the  exception  of  Just  tvo  of 
the  i  inp:tB,  It  'c  aasuiaed  here  that  tlxe  ciuantitles  of  rsw  gasolines  avail¬ 
able  ruul  thf-t  the  quality  clamcter isllcs  of  tliese  materials  are  Independent 
of  decisions  taken  in  the  blsixling  department.  In  fact,  from  the  vlevpoint 
of  a  i-efinery  superintendent,  the  r^uantlties  and  characteristics  of  all  the 
Gti-earu3  arc  variables  subject  to  ccintrol.  By  altering  the  choice  of  crudes, 
t.he  I'eactor  teuperaturea ,  the  recycle  ratios,  and  the  assignment  of  Inter- 
uedlatf'  oil  streamc,  the  central  managerwnt  is  able  to  Influence  the  slxe 
ajKl  conrpcc': tlon  cf  tl*o  gacolaie  blending  materials.  This  type  of  Interlock 
can  I  iHy  be  studied  in  a  lai-ger  uvAiel  that  cuts  across  departmental  lines. 

A  I'cc  T*.i^pci;;t  is  the  assunxptlou  that 

adclUi  nsJ.  .quantities  of  finished  refinery  products  may  bo  marloBted  at  the 
prevailir.g  reallsat ! ons .  In  t/.e  casr  jf  a  major  refiner,  this  is  open  to 
serious  dcubt.  It  is  n..'t  at  all  obvious  trjit  a  large  cciupeny  can  put  additional 

quantities  on  the  mBiliPt  vlt'.iout  disturbing  the  current  price  structure.^ 

7 - 

During  the  course  of  short  tine  intervale,  the  quantities  and  quaillties  of 
the  available  blend  IrtT  stocks  any  change  radically.  An  unforeecen  shutdenm 
of  a  unit  can  ta>:e  place.  A  new  unit  nay  be  breu^t  on  atream.  A  n«v  source 
of  crude  oil  may  be  introduce!,  etc.  Fr>r  these  reasons,  in  order  to  keep 
blending  ansdyses  up  to  date,  it  is  the  current  practice  vlthin  the  company 
to  recalculate  schedules  nt  leant  ^nce  a  month.  The  ptiroia:  aliocntiona  cf 
raw  streams  among  the  various  uses  could  be  significantly  affected  if,  say, 
a  fire  occurred  in  a  large  catalytic  cracker. 

See  J.  S.  Bain,  The  BcododIcb  of  Uie  Pad  fic  Ooaat  Petrol  cun  Industry, 
University  of  Oalifomia  Presa,  Part  t,  l  Al,  Part  II,  19^5;  and  Part  III ,  1^7. 
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At  the  B&me  tlnve,  :io  one  —  least  oT  all  an  oconondst  —  Is  In  a  position 

to  make  a  reasonable  quantitative  estimate  of  the  relationship  between  a 

company’s  price  structure  anJ  the  demand  for  its  products. 

For  purposes  of  computation,  an  explicit  marketing  limitation  could  have 

been  placed  upon  each  finished  pro-luct.  In  onier  to  provide  a  benchmark, 

and  then  one 

however,  it  was  decided  to  study  one  raslc  case/^ln  which  a  sales  restriction 
was  introduced.  The  comparison  of  these  two  provides  an  indication  of  the  ex¬ 
tent  to  which  profits  could  be  Improved  by  overcominp;  the  marketln;»  harrier. 

The  comparison  also  Indicates  the  ma/rnltude  of  the  I’ap  between  the  reflrier;/ 
realisation  and  the  incremental  cost  In  tne  second  Instance, 

Correspondin'’  to  the  assumption  about  the  realisation  upon  refiner^'  pro¬ 
ducts,  it  was  ai^roed  that  one  of  the  inputs  —  tetraethyl  leal  —  would  be 
considered  as  available  without  any  quantity  restrictions  at  the  Ethyl  Corpora¬ 
tion’s  prevailing?  tanx-car  lellvorea  price.  During?  wartime,  TEL  bas  l^^en 

rationed  amonr  refiners,  but  at  least  during  1953 »  this  type  of  rnstriction 
hold  true.  the 

did  not4  If  lead  were  to  become  short  a,’aln,<^ 

explicit 

availability  limitation  would  have  to  oe  considered  hv  anAelement  in  the  proole 
As  shown  in  the  prece’.ln'?  chapter,  it  is  possible  for  the  refiner  to  substitute 
reformed  and  cracked  <?asolines  over  a  wide  ran  ;e  in  place  of  tetraethyl  lead  in 
the  finished  prcxlucts.  Although  lead  is  employed  in  virtually  every  gallon  of 
ganoline  today,  its  use  in  motor  fuel  could  be  ellmlnate'd  —  but  only  at  a 
distinct  cost  in  dollars  and  cent#  to  the  refiner. 

Havin'?  decided  upon  the  two  ma.ior  assumptions  —  fixed  ^^ssoline  stock  avail 
abilities  and  fixed  refinery  product  realizations  —  the  next  step  was  to  de¬ 
termine  which  of  the  numerous  quality  characteristics  of  the  finished  products 


would  b«  studied  explicitly.  The  MsoufacturiJig  Scoacxaics  Dl>ri8ioxi  pointed 
out  that  there  were  oany  specificatians  that  finished  laotor  fuel  sust  meet  -- 
three  types  of  octane  numbers;  lO^t,  3t4,  and  end  points  of  the 

boiling  range;  sulfur  ajod  gua  content;  vapor  pressure,  etc.  Any  of  these 
could  conceivably  limit  the  output  of  either  prealum  or  regular  grade  gaaollne 
For  purposes  of  this  study,  the  ones  that  seeaed  to  be  worth  the  most 

close  attention  were:  Reseeirch  octane  nuaber,  the  percent  sulfUr  content,  and 
a  volatility  Index.  In  addition,  the  end  point  problem  was  handled  auto¬ 
mat  icmlly  by  limiting  the  end  points  of  the  individual  raw  gasolines.  During 
the  course  of  the  counting  runs,  it  turned  out  that  eulfur  was  not  a  limiting 
specification,  but  that  the  octane  number  and  volatility  index  specifications 
did  OGM  into  pLaj. 

In  order  to  protect  proprietary  inforaatioa,  the  emct  hoost-brtnd  sped- 

flcatlTns  cannot  be  revealed  here.  II  can,  Licwever,  be  said  that  the  r6CX) 
premium  gasoline  has  a  hlgJier  nlnlinum  octane  number  and  volatility  index  -- 
and  a  lower  maxi  in nn  pennleslble  sulfur  content  and  end  point  than  the  regular 
76  gasoline.  On  both  products,  the  mojclr'ua  pennlBSlble  tetraethyl  fluid 
concentration  la  3-0  olll lliters  per  gallon. 

t 

Although  the  blending  section  rust  operate  within  whatever  cenpany  speci¬ 
fications  arc  currently  effective,  the  reader  should  not  Bujipcse  that  these 
limits  remain  fijeed  over  time.  At  any  noaent,  "octane  wars"  between  refiners 
can  break  out.  Ihc  effect  upon  motorists  --  and  upon  refinery  profits  --  is 
of  much  the  saxae  nature  as  an  actual  price  war.  For  this  reason,  a  separate 
coraputatloo  waa  set  uj.  in  order  to  determine  the  effects  cf  a  2.5  octane 
nuid>er  increase  in  the  specif! cat  ions  of  tl.c  preLn'-um  gasolln'’. 
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Before  proceeding  to  the  technolo#^  of  the  problem,  two  more  of  the 

economic  eseumptlons  should  be  brought  out  —  one  hevln/?  to  io  vlth  the 

are  manufactured, 

cutter  stocks  and  the  other  wi  th  transportation  costs ,  The  raw  gasolines 

and  the  blending  operations  are  performed  at  two  distinct  locations  -- 

at  Wilmlngt  on  In  the  Lrs  Angeles  area  and  at  Oleum  on  San  Francisc)  Hay. 

Raw  stocks  are  continually  lifted  by  tanker  for  interchange  between  the  two 

refineries.  Union  Oil  estimates  its  tanxer  costs  between  the  two  points 

to  be  nominal  -  only  a  few  percent  of  the  refinery  realiaation  on  prcxi’am 

gasoline.  In  the  analysis  of  the  coinpar\y*s  gasoline  blending  problems, 

it  is  evidently  a  desirable  thing  to  take  account  of  transportation  costs 

and  of  differences  in  product  realisations  at  the  two  refineries.  At  the 

same  time.  It  was  ^^aay  to  see  that  by  bring’ ng  geography  into  thf  problem, 

the  volume  of  computations  would  have  increaael  by  more  than  100^.  The 

storage  capacity  of  the  then  available  computing  laachlne  —  the  International 

business  Kachlne  Card  Progranraed  Calculator  —  would  have  been  completely 

Inadequate  to  take  care  cf  the  enlarged  model. As  a  first  effort,  therefore. 

It  was  agreed  that  the  geography  should  be  left  out  of  account, 

conaldered  to  be 

raw  gasoline  stocks  wercAarallable  in  one  company-wide  pool, 
were 

and  the  finished  pro'iucte/^  valued  at  the  compjany’s  overall 

realisations,  (Roughly  speaking,  these  revenue  flguren  are  of  the  same  order 
of  magnitude  as  the  current  spot  market  refinery  tank-car  quotatlonn. ) 

4*  This  type  of  limitation  would  not  have  existed  If  one  of  the  higher- 
capacity  electronic  computers  had  been  available,  e.g,,  the  IBM  7bl  or 
the  Univac, 


i 
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th«re  Is 

fljMLljr,  tbm  problsa  of  ai2tt«r  stock  rUlMtfoai.  As 

Indicated  by  Flfurs  10  ,  ei^it  of  tbs  hesnry  psol  inw  Biigr  be  tilflQfifl,  not 
oeiy  for  motor  ftiel  oon^oeeiitib  but  also  fbr  blandlnt  vlth  hscvy  rsslduai 
into  lo.  6  fuel  oil.  Ihe  flfure  also  shove  that  this  besTy  rssiduua 
be  eomrorted  Into  lo.  6  oil  by  the  sltenMtlvs  route  of  blending  with  So.  2 
fVtmaos  oll.^  Qtader  the  sssis^itlon  that  set  heavy  residuua  Is  alsiqns 
going  to  be  cot  lava  to  the  visoosity  qtadfl  oat  loos  by  blending  vlth  So.  A 
oil,  it  is  possible  to  establish  a  valw  for  the  heavy  resldut,  sad  tht  to 
go  oa  to  datenlns  a  cutter  value  for  any  particular  heavy  gasoLit. 

To  give  a  yuroly  hypothetical  illustmtioD  of  the  prooedurs,  asatavt 
the  followlag  product  valuesi 


So.  2  funaoe  oU  $h.OO/^. 


So.  6  fual  oU  $2.00/b 

^  of  lo.  2  oil  in  Ho.  6  blend  reared  for  cutting  bs 
to  lo.  6  oil  speelfleatlans 


heavy 

sr 


$  of  heavy  gaadine  la  lo.  6  blend  required  for  oattlag  bach 
heavy  realduua  to  So.  6  oil  apeelfloatlene  1^ 

If  No.  2  oil  is  used  for  branding,  the  valuo  of  the  heavy  rosiduuB,  p^, 

may  be  established  through  the  fo^i.airliig  equations 


faint  at  Ho.  6  oil  -  falue  of  Ho.  2  ♦  faius  of  raaHiiT 

oU  oosvottBBt  nr^iri— ♦. 

$2.00  -  .25($»i.OO)  ♦  (1  -  .25)(Py) 

p,  -  $a..33y$ 

gsaarsi,  the  cutter  vUl  not  he  a  ftnIShed  So.  2  oil,  but  ratbar  eoas  gaa 
oil  ossveaaat  that  aoald.  ha  uead  for  lo.  2.  Ihas  these  are  Hvwal  altaraa 
tlve  gaa  oils  that  wn  ka  sig^liiged  for  euttar,  the  prOhInsi  tahsa  mi  adilti<n»l 
rmsplantty.  Sm  Chapter  fl  belov,  p.  VI  -  26 


L  A  ond 
0(«ufrt 
products 


Mo'xetoble 

products 
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Next,  this  velue  for  is  substituted  in  e  costing  equation  for  the 

heavy  gasoline  cutter  stock.  Out  of  this  relation,  there  is  derived  P  , 

the  cutter  value  of  the  particular  gasoline. 

Value  of  No.  6  oil  »  Value  of  gasoline  ^  Value  of  residuum 

component  component 

12.00  .  .20(p^)  r  (1  -  .20)(J1.333) 

p  *  %U.668/B 

D 

Values  for  each  of  the  eight  possible  gasoline  materials  were  calculated 

in  this  manner.  Certain  treatment  cost  savings  also  had  to  be  credited  to 

the  cutter  operation.  The  final  calculated  net  cutter  values  range  from 

6P.2f  to  P7.^  of  the  refinery  realisation  on  premium  gasoline.  In  effect, 

this  is  an  opportunity  cost  calculation  of  the  same  nature  as  the  ** gasoline 

replacement”  technique  described  above,  pp.  II-5-7*  Like  the  gasoline 

replacement  fonsula,  this  method  applies  to  the  comparison  of  Just  two 

alternative  refining  operations.  (In  the  one  case,  the  alternatives  are 

fuel 

cracking  versus  fuel  oil,  and  in  the  other,  gasoline  versus  No.  2Aoil 

were  capable  of  meeting  the 

cutter  materials.)  In  the  event  that  the  company’s  residuur/y 
fuel  oil  specifications  without  any  additional  cutter,  the  heavy  gasoline  could 
not  oe  given  the  premiuir  value  as  cutter,  but  would  have  to  be  assigned  a 


lower  credit 


2.  The  blending,  tachnolo/y 
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Of  the  three  relevant  (gasoline  properties,  two  —  the  volatility  index 

and  the  percent  sulfur  content  — -  rise  to  straifjht forward  linear  blendlnp; 

relationships.  For  example,  in  a  blend  of  two  raw  j^asollnes  with  sulfur 

contents  of  O.lCrt  and  0,yiif  the  resulting  mixture  will  contain  0.20'f  sulfur. 

Similarly  with  the  volatility  index,  a  blend  will  have  an  index  that 

lies  midway  between  those  of  its  two  component  raw  '’:asollne3. 

Unfortunately,  the  octane  nuriber  of  blends  cannot  be  prelic"^ed  in  such 

a  direct  manner  from  the  octane  numbers  of  the  components.  Like  most  kinds 

<ir. 

of  enjilneering  work,  the  forecast  of  .gasoline  octane  numbers  is/v 
rather 

artAthar,  an  exact  science.  Refiners  themselves  necessarily  make  pape** 
predictions 

^of  the  octane  number  of  proposed  blend:?,  but  before  marketing  a  product 
they  will  alir.ost  Inevitably  taxe  the  precaution  of  testing  the  m.Lxture  in  an 
octane  ratlnr  engine. 

One  phase  of  the  problem  has  he<»n  investigate  i  extensively  by  the 

refining  industry  —  the  relationship  between  octane  numbers  and  the  TEL 

concont ration  of  a  particular  .-asoline.  Just  as  in  the  rapfitha  reforming 

problem  above,  it  has  been  observed  that,  as  the  lead  concentration  level 

is  increases,  the  octane  number  increases  at  a  decreasing  rate.  In  order 

to  predict  this  relationship,  there  is  iri  widespread  use  tcnday  %n  othvl 

blending  nvnograsi  put  fo"warj  by  Hebi,  Reniel,  and  Gartor  in 
fortunate '.y. 

There  is,  Na  throe-parameter  analj-tlc  expression  that  gives  a 

close  appro. irrat lor  to  the  results  predicted  by  this  blending  chart,  ^ere  t 


.  L.  E.  Hebl,  T.  B.  Renriel,  ami  F.  L.  Garton,  "Ethyl  ‘•’luid  Blemiing  Chart  for 
Motor-Methmi  Octane  Numbers",  Industrial  and  Engineering  Chemistry.  Vol.  31, 
July,  1939,  pp.  %2-b^5. 
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representa  the  octfLne  nunl>er  of  the  leaded  jjaaollne,  and  x 


TBL  concentration  level  (in  alllllltcrs  per  gallon): 


(Y.l)  t  ^  a  +  bx 


the 


For  aay  \cikxuam  gaaoline,  the  U.rce  positive  paramators,  a,  b,  and  c, 

have  to  be  drten&iJied.  any  be  accoopllshad  readily  by  taking  three 

the  reau.tin^’ 

obMrvaticna  of  octane  nuinber  for  various  lead  levala,  and  solving  /\ 

aljaultaneous  linaar  equatiooa.  Using  equation  (Y.l)  in  this  way,  sets  of 

calculated 

eonatanta  verey^  for  all  19  of  the  Ikilon  Oil  Coopany'c  rav  gasolines. 

Testing  aevan  lead  levels  in  each  of  these  I9  cases,  the  equation  has  virtuailly 
aleaya  given  a  prediction  that  lias  vlthin  -  0.2  octane  msaber  of  the  one 
yleldad  by  the  blending  chart. 

No  specisd  significance  should  be  attached  to  the  particular  fom  of 
equation  (V.I).  It  la  aerely  a  device  for  enabling  a  ccc^xjtlng  aachine  to 
parfom  the  aaae  csLlculation  as  a  refiner  vrith  his  ethyl  chart  and  his  straight 

edge. 

A  secGod  phase  of  the  blending  problera  is  nore  controversial  than  the  TEL 
aspect.  For  a  lead  Level  of,  say.  3'U  nl.,  and  for  a  90-50  mixture  of  tvo 
gasolines  A  and  B,  refiners  frequently  calculate  the  cetane  nuaber  of  the  blend 
to  be  the  30-90  vslghted  average  of  the  octane  nuaber  of  gasolines  A  and  B  -- 
each  with  3-0  al.  l%e  weighted  average  iqijpears  satisfactory  for  oan^  stocks  -- 


especially  whan  both  ccaqpocents  are  of  a  psuraffinlc  nature.  There  arc,  hcvever, 

7 

at  least  tvo  papers  puhlicly  available  --  one  by  Baatoan  and  the  other  by 


Bogen  and  Nichols  --  that  call  into  question  the  otralght-line  average  method. 

7  - 

Du  Bois  Eastoan,  "Prediction  of  Octane  Nuatbers  bdA  Lead  dusccptibllities  of 

□siaoliDe  Blands",  Industrial  and  Engineering  Oaealrtry,  Vol.  33,  Dccei:±>er,  l^^l, 
pp.  1553-1560. 

8 

J.  S.  Bogen  and  R.  M.  Nichols,  Oalculatlng  the  Perfonaance  of  Motor  Fuel 
Blends",  Industrial  and  lagineering  Chemistry,  Vol.  hi,  November,  19^  S  PP- 

2629-2635. 


;  ;  ^ 

:  -A*' 

Both  papers  iDdlcate  that  as  the  percents^  of  the  high  octane  coBq>onent  In 
a  binary  lalxture  Increases  (both  coo^xanenta  initially  at  the  aame  TSL  con- 
oentration  level),  the  octsAe  level  of  the  nixture  may  increase  at  a  decreasing 
rate.  In  other  words,  the  octane  number  of  the  blend  tends  to  exceed  the 
weighted  average  of  the  octane  nuraber  of  the  two  conponenta. 

Bogen  and  Ifichols  are  prliarlly  concerned  with  three -coraponent  mixtures, 
and  they  report  only  five  observations  at  each  of  three  Ttt.  levels  for  strictly 
binary  blends.  The  set  of  15  points  has  been  plotted,  and  Is  reproduced  as 
Figure  II.  For  each  TEL  concentration  level  shown  on  this  chart,  strai^t- 
line  interpolation  hsis  been  used  between  adjacent  obaervatlons.  In  all  like- 
lUujod,  these  straight  lines  also  underestimate  the  octane  number  that  is 
attainable  by  blending.  The  existing  data,  unfortunately,  do  not  provide  a 
basis  for  curTlllnaar  Interpclation. 

At  first  glSLOoe,  it  would  appear  to  be  a  trivial  aatter  whether  octane 
nuo^rs  were  calculated  by  drawing  straight  lines  between  two  adjacent  points 
or  between  two  end  points  of  e  constant-TIL  line,  nevertheless,  a  eaiyile 

Q 

calculation  presented  in  an  earlier  paper  indicates  that  the  two  methods 
can  lead  to  T*L  input  requirements  that  differ  by  20^. 

It  Is  evidently  worthwhile  for  a  refiner  to  incur  considerable  expense 

in  testing  gasoline  blends,  and  in  attas^lng  tc  predict  the  occurrence  of 

curvilinear  octane  relationshlj s .  Nevertheless,  at  the  present  stage,  the 

Union  Oil  Company  did  not  have  such  blending  data  available,  and  it  was 

necessary  to  revert  back  to  the  more  conventional  straight-line  interpolation 

procedure  for  calculating  the  octane  numbers  of  blends. The  information 

^ .  ■  ■  ■  '  ^^})holo-of fsst 

A.  8.  Hamne.  "Concave  Prograssning  for  Qasoline  Blrfnds",  P-393 

paper),  The  RAHD  OorporatioD,  Bsurta  Nosiioa,  California,  April  23,  1953- 

1C 

In  going  back  to  the  sisgle  straight-line  octane  interpolation  procedure, 
this  paper  follows  ths  model  proposed  by  A.  (Thames,  W.  w.  Cooper,  and  B.  Melljn, 
"Blending  Aviation  Qasolines  --  A  9tudy  in  Programing  Istsrdspexwlent  Activities 
in  an  Integrated  Oil  Cos^mmy",  Iccnometrica,  Vol.  20,  No.  ?,  April,  1932,  pp. 
133-159.  But  unlike  the  avlatlao  gasolinTaysten,  the  preasnt  one  includes 
ethyl  fluid  concentration  levels  as  a  specific  variable  with  a  curvilinear  effect. 


%  fher.TQii^  cracktri  gasoH'*? 


t '  1 0  f  y  blend 


Fig. I  i -Octane  number  versus  percent  of  therrpally 
crocked  component  in  binory  blend.  Three  tEL-l€fvels 


m 
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curr«ntly  at  hand  made  it  possible  for  the  analysis  to  take  account  of 
non-Unear  TEL  relationships,  but  not  of  the  non-linear  p^s^Hne  blendlni» 
features . 

The  TEL  inputs  are  one  of  the  elements  in  the  Union  Oil  Company 
problem  that  involved  a  departure  from  a  strictly  linear  system.  In 
addition,  the  cutting:  temperatures  on  two  of  the  /gasoline  distillation 
units  introduced  non-linear  effects.  For  purposes  of  blending,  just  as 
in  the  Kelloc^  Company  naphtha  refomiin>:  problem,  the  cut  point  on  each 
of  these  units  is  an  independent  variable  that  is  subject  to  control. 

The  cutting  temperature  influences,  not  only  the  absolute  qiiantity  of 
material  that  is  available  for  motor  fuel  blending,  but  also  the  quality 
of  this  component.  Figures  12a  -  12c  contain  the  necessary  data  for 
one  of  these  two  distillation  units  ~  the  one  producing  gasoline  stock  #31» 

According  to  Figure  12a,  as  the  cutting  temperature  moves  from  300® 
to  /.OO®  F.,  the  amount  of  material  available  for  motor  fuel  blending  rises 
steadily  from  1,900  3/CD  to  3,050  B/CD.  (For  each  additional  barrel  of 
motor  fuel  obtained,  there  is  approximately  a  one-barrel  decrease  in  the 
amount  of  gasoline  cutter  stock  available.) 

curves, 

Charts  12. b  and  12. c  contain  the  original  Union  Oil  Company  A 

a^so  indicate  the  fom  of 

and/^the  equations  that  were  fitted  to  these  lata.  As  the  make  of  motor 

virtually 

fuel  increases.  Figure  12. b  indicates  a/|linear  decline  in  the  octane  number  — 
both  clear  and  with  3  nil.  of  TEL/guHon.  Also,  as  the  make  of  motor  fuel 
increases,  there  is  a  steady  linear  oecrease  in  the  volatility  index,  and  a 
quadratically  increasing  sulfur  content.^ 


11.  AlthoUfii  when  distilling  this  particular  stock  the  sulfur  content  increases 
quadratically,  on  the  other  material  (gasoline  #36)  the  sulfur  content  only 
changes  linearly. 


■9  -vyr*. 

V’-i5 
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■Clear  and  leaded  octane  numbers  vs  x^, 
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3.  The  ■athenmttcal  ■odel 
3.1  The  basic  form 

lOuhn  and  lUeker  have  proved  the  follovlog  "equlvaleace  theorem": 

"Let  the  functloxie  f  (*)»  g(*)  ^  concave 

aa  well  as  differentiable  for  x  >0.  Then  Is  a 
solutioQ  of  the  nBixinuin  problem  If,  and  only  if,  x° 
and  saae  give  a  solution  of  the  saddle  value  problem 

12,  13 

for  0(x,u)  s  g(x)  u'rx." 

In  other  words,  they  have  established  that  the  solution  to  a 
concave  constrained '^axloua  problea  is  equivalent  to  a  certain  miui- 
imx  aaddlepoint .  At  the  suggestion  of  Harry  Markowlti,  tlie 

gasoline  blending  aodel  was  set  up  in  the  appropriate  minaax  form. 

An  Iterative  digital  technique  was  then  devised  for  solving  the  problem.  ^ 

Fron  a  rigorous  ntheoatloal  standpoint,  the  present  coegHitlng  procedure 
is  an  unsatisfactory  one,  but  the 

^  method  does  lead 

to  solutions  that  should  be  sufficiently  accurate  froa  the  standpoint 
of  refinery  sanagere.  Convergence  tabes  place  without  an  inordinate 
aaount  of  computing  time.  All  calculations  were  performed  on  a  machine 

The  expression  u'Fx  is  shorthand  for  the  folloving: 

Uj^  fj^(x)  ♦  Ug  fjCx)  ♦  .  .  .  ♦  u^f^(x).  The  are  non-negative 

Lagrangsan  ailtlpliars,  and  correspond  cloaaly  to  the  econoadat's 

notion  of  "shadow  prices". 

V.  Kuhn  and  A.  V.  Tucker,  "Monlloear  Programing",  Proceeding  of 
tte  fcc<^  Berlaty  >Byaiua  on  Msthaaatical  Statistics  end  ^■ohaED.lTy, 

517  J*.  l^yiM,  dilvsnity  of  Chlii^omia  Wmss,  p.  535. 

U 

f^kowlts  has  experiaentad  with  an  electric  analogue  solution  for  this 
class  of  applications.  The  present  blending  aodel  involvlnc  19  gasoline 
stocks  would  clearly  svaa^  the  capacity  of  any  of  the  present  generation 
of  aarhtnas. 
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thmt  is  readily  available  throuf'hout  the  United  States  —  an  IBP 
Card  Prograncned  Calculator,  Yodel  2. 

Dr.  Martin  Beckman  has  objects  i  to  the  application  of  the  Kuhn- 
Tucker  theorem  in  this  particular  instance.  He  has  made  the  point 
that  for  a  function  of  raary  variables,  f(x),  to  be  truly  concave.  It 
is  necessary  but  not  sufficient  for  that  function  to  be  concave  in 
each  variable  taken  separately.  It  has  been  possible  to  devise  ex¬ 
amples  that  satisfy  all  conditions  of  the  .-asoline  blending  probler,  — 
that  is,  concave  in  each  variable  taken  separately  —  but  which  never¬ 
theless  violate  the  strict  Kuhn-Tucker  con'lltions.^^  The  cour.ter- 

15.  The  following’  example  is  not  derived  from  any  Union  Oil  Company  data, 
and  serves  only  for  illustrative  purposes.  Consider  two  blending’ 
stocks,  labelled  "A"  and  **H"  with  the  following  characteristics: 


Stock  A 

Stock  B 

Octane 

number,  clear 

74.0 

90.0 

Octane 

number,  with  3»0  ml. 

of  TEL  90.0 

Number  of  lal./gal.  required 
*?C  octane  leader  gasoline 

for 

3.0 

0 

Assuming  straight  line  interpolation  at  constant  lead  levels  for  the 
two  stocks,  and  employin’  the  Hebl-Rende 1-Carton  blending  chart,  it 
will  be  found  that  a  50-50  mixture  would  require  1.6  —  and  not  1.5  — 
ml. /gal.  in  order  to  reach  the  90.0  octane  level.  This  clearly  consti¬ 
tutes  a  violation  of  the  concavity  assumption,  but  a  refiner  will  also 
recognize  that  the  example  requires  an  unusual  configuration  of  octane 
numbers  and  lead  susceptibilities. 


h 


exaaples  glre  rise  to  nultlple  Isolated  ■axiisa.  c  ncavlty 

it  is  possible 

property  in  the  vmrlouB  functions,  ^tr  adopt  bleixiing 

•chedulea  that  are  locally  optliaal,  but  which  do  not.  In  fact,  repro- 
aent  the  beat  of  all  aolutions  available . 

Ibe  problea  cannot  be  dlsmlased  lightly,  but  I  suspect  that  in 
■oat  partlCTilar  applications  sucl.  local  extrcaa  can  be  detected  by 
the  Inveatl^tor.  The  difficulties  arising  'ut  of  n  ?n-cwncnv:  *  y  ar'.* 
not  unique  to  the  progrannlng  method  describ*^!  here.  The  saac  argicaent 
would  apply  If  the  TEL  levels  were  taJcen  as  parameters,  and  a  ^;»•i■«l^ate 
lineaur  programaing  problea  tlien  run  cdY  for  each  of  raan,>'  coablnatl.  .ns 
of  the  two.  Thla  difficulty  would  cone  uj'  in  the  case  of  aiy  trLol- 
aad -error  aathod  that  conalsted  In  changing  one  variable  at  a  time, 
and  observing  the  effects  upon  the  payoff.  And  finally,  ty.e  failure 
of  the  proper  curvature  conditions  could  frustrate  optlnizatlons  by 
way  of  a  I^mge-Lemer  ahadov  price  laarkct  mcchaninn. 

3.2  Definition  of  independent  ^/arlablee 

It  now  reaains  to  set  up  the  fonaal  relatlonahips  cc  rrespondix^ 
to  the  FDohn-Tucker  payoff  function,  g(x),  and  t];c  set  of  restraints, 
fj^(x).  Table  V.l  presents  achcaatlcaliy  the  Independent  non-ne^- 
tlve  Tarlahlas,  Xj  ,  which  mist  be  selectci  in  .  rdcr  to  Jet-':T:ir.r  th.r 
bLendlag  schedule.  The  whole  set  of  23  will  be  termed  the  vector  (x). 

There  are  1  ■  raw  gasoline  stocks,  labelled  3I.  3'.  ,  bl,  ...  80 ,  as 
indicated  by  the  row  headings.  In  each  Instance,  we  are  given  ,  the 
total  quantity  of  this  natcrial  available  for  blending.  Ti.c  problem 
■naista  of  assigning  various  aamunta  of  the  raw  gasolines  ano  ng  the 
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Tikble  V.l  OMoline  Blending  Problea  --  Definition  of 
Various  Quantities 


Stocks 

1 

_ 1 

Quantity  of  stock 

J  assigned  to 

Udentix'i- 
cation  Bo. 

^lantities  | 

Available  I 

thiw  pirxluc 

ts,  B/CD 

b/cd  i 

(1)  7600  pi'enium  (2)  76  regular  (3)  Cutter 
motor  fuel  motor  fuel  Stock 

31 

*^31 

3,050  ! 

0 

*31 

^'^"*31^ 

36 

^36 

2,060  ! 

j 

0 

*36 

''‘36-V’ 

Ul 

1,000  • 
i 

*41 

0 

U2 

2,910  * 

1 

^5 

*42 

^S*2**42"*45^ 

43 

^3 

3,460  \ 

1 

*46 

*43 

K3’*43"*46^ 

44 

'»44 

1,730  : 

i 

0 

*44 

^S^"*44^ 

51 

^*51 

no  1 

*51 

0 

52 

S2 

1,240  i 

*55 

*52 

(S2-S2-^5> 

53 

*^3 

860  I 

1 

*56 

*53 

^‘*53‘*53‘*56^ 

54 

'»54 

490  1 

1 

0 

*54 

tw> 

60 

Si) 

6,150  1 

^^■*60^ 

*60 

0 

6I 

3,1*30  1 

^Si"*6i^ 

*61 

0 

62 

1,460  1 

^  S]2"*62^ 

*62 

0 

63 

3,200 

*63 

0 

64 

1,600  ( 

^‘*64**64^ 

*64 

0 

65 

j 

n,230  j 

^‘*65**65^ 

*65 

0 

66 

%6 

r,440  i 

<'>66-»66> 

*66 

0 

70 

**70 

1,700  j 

j 

^‘*'ro"*70^  *  **70 

*ro  ■ 

0  0 

do 

0,710  ! 

] 

'^eo-W  =  ° 

*80* 

’60  ° 

Totals 

— 

61,910 

i 

*1 

‘2 

— 

1 

TIL  concentration, 
(■l./sal.) 

*21 

*22 

three  alternative  viscb  --  (!'  '’'Vv",  ’'rcnl-'j-!  grade  nmtor  fuel,  (?)  '6, 
reg'uier  gntdfc  aot'^r  fjel ,  -rd  (3)  cutter  stock.  In  addition,  Lt  l6 


ncccGsnry  to  oeleot  x  .  and  x,„^,  respect  Iv^'ly ,  the  ethyl  Pnilti  enn- 

2J. 

centratlons  In  preralun  and  regular  gasoline. 

The  ccaaixitatlona  ware  not  aet  up  with  the  full  generadlty  of 

considering  I'’ 19)  ''V  i-udlvldual  possible  gasoline  components. 

Inctoed,  for  'Tjr''oecs  of  c  :or.orilzing  on  nech'r.o  oapacitv,  it  whp 
to 

agreed^  rake  a'A’an^c  gueQ3''o  at  vdii'^h  of  the  elenente 

would  turn  ort  to  have  zero  •mlvies.  (Wherever  a  component  was  pre- 

aaclgned  In  this  way.  tliere  Is  a  zero  entry  In  the  corresponding  box 

of  Table  V.l.)  It  was  also  Bgree<l  that  It  would  not  be  econcmlcal 

16 


to  thro/  away  any  of  these  particular  gasoline  atockc.  In  this 
manner,  It  was  ]K)iBtblc  to  reduce  the  number  of  Independent  gasoline 


17 

cooiponento  from  ‘  '  down  t'^  21. 

The  Lagrongean  rmdtlplloro ,  ,  taken  directly  frum  the  basic 

cal culat lens .  enable]  us  to  clieck  up  on  these  pair ticular  aaouaptlons 
at  a  later  stage  in  the  analysis.  Fortunatcl:i^ ,  fne  initial  guesses 
turned  out  to  be  correct  or^es.  Despite  the  good  luck  for  the  pn.rt.1  mlar 
seta  of  poraaeters  studied,  these  aasuarptlons  -were  an  'insat  Is  factory 
short-cut,  and  can  be  dust  1  fled  only  on  grounds  of  connjtatlnnaLl  con¬ 
venience  . 

15 - 

Sven  If  a  material  could  not  be  employed  directly  as  cutter  or  as 

regular  grailc  gasoline,  under  the  1953  price  atruciure,  It  would 
always  pay  to  degrade  book  flnisheu  ('tTOC  gasoline  for  purposes  of 
blending  the  ^urtlcular  material  up  to  the  specifications  for  regular 
grade  gasoline . 

17 

For  syasactry  In  the  ‘^ubf;^';uent  work,  there  are  t-vc  additional  constants 
that  are  written  as  If  they  were  Independent  variable*.  Tdic  tem  x.^ 

Is  Identically  zero,  and  Identically  equals  Cj^  . 


Table  V.I  alwiiyB  refers  to  the  quantity  of  stock  J  in  regular 
gr«*lo  gasolloe  as  xj.  Even  In  the  caae  of  stocks  ^3^  1^*3'  >  where 

the  cutting  teraperature  can  equally  well  be  regarde<l  as  the  Independent 
variable,  this  saae  convention  was  Cinployed.  (Cu  Figuro  12. a.  note 
the  one-to-one  relatlcmahlp  between  the  cut  point  and  the  quantity  of 
notjr  fuel  avaliabl''.) 

To  slapllfy  further  referent,  tv't  new  dependant  variables  are 
also  deflnrsi: 

(V.2)  q  -X  )  +  ,t  ■=  total  production  rate, 

^  B/CD,  for  preaivua  gratie 

(^aol  Inc 

(J  X  Ul,  ro-ou,  ;o) 

(k  ^  Uv.  fvo.  -s  /) 

,V.-)  T.^  X.  ^  total  production  rate, 

J  b/CD,  for  regular  f^rade 

gaaol  Lne 

(j  31.  ,  -’iiJl.  M-',;*.  o-t^. ,  d^-) 


3  •  j  Dcflnltijn  of  the  payoff  rmctlon  and  ivstralnts 
The  neyoff,  6(x).  ir.  determined  by  crediting  the  gross  refinery 
reallzatl’in  fron  the  ael^  of  pre.-alun  and  regular  grsule  gaaollnc,  and 
then  charging  up  t'r^  ctlv'l  fluid  coste  and  the  loss  of  7>otentlal  cutter 
credltf  acBoelate^l  with  the  partlculAi'  blcndlnq-:  schedule: 

(v.'O  «(X)  ,  (,.^-r^,  *  rjTp,, 


P:,?  t' 

^  -  X 


J  ’  T1 


o 


In  equation  (V.U),  the  unit  costs  oikl  prices  are  iei'lned 


8.6  follovs: 


(V.‘^)  refinery  real! ratlin  '■n  T'OO,  premivun 

gBLfloline,  $ 

*  refinery  real  1  nation  on  T^,  regi.ila.r  gra^lc 
gaaol  in' ,  $  /B 

-•  cost  of  ethp'l  fluid,  $/U2  ml. 

Pj  unit  cutter  cred.lt  on  st-.ck  J, 

;j  -  31.  36,  h2,  U3,  U,  S2,  ^3,  ^^0 


/ 


The  blending  schedule  x)  nsi£t  not  call  for  the  pn>iuction  of 
negative  amounts  of  any  motor  fuel  or  cutter  ccnponent ,  ajxi  so  there 
are  certain  upper  aivi  lower  Units  Imposed  upon  the  various  Xj.  In 

addition,  the  schedule  must  not  require  TKL  concentration  levels 

exceo  Jinj, 


3.00  ml. /gallon.  And  because  of  the  flash  point  specif  1  cat  1  ui. 
on  finished  No.  (1  oil,  It  must  not  require  cutter  conponents  #31  or 
#30  to  fall  below  q^  and  q^^^ ,  respectively. 


0  • 

*21  - 

3.G 

0  ■ 

*22  - 

3.0 

J 

■^u  - 

- 

^31 

0 

4<,  - 

s.  - 

0 

i. 

H 

(1  -  41.  1*4,  11,  :!*, 

( J  -  1*2,  43.  4‘  .  4t,  S2,  13.  'o,,  ) 
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For  stocks  U3,  52,  and  53,  the  quactities  of  individual 

cutter  stocks,  laist  not  be  negative; 

(V.7)  f^(«)  .  >  0 

V*)  ■  Hi  -  =S.3  -  2  0 

r^(x)  =  -  X5J  .  .55  >  u 

»  1,3  -  -  Xjg  >  0 

Tba  six  I’SBilnlm  •<iaxtloiis  are  eoBiisetad  vlth  the  three  epeeifloe* 

tlons  on  each  of  the  two  motor  fuel  products.  The  functions  f^(x)  and 
fg(x)  relate,  pcspectlvely,  to  the  octane  nuober  of  preTilun  and  regular 
grade  taotor  fuel;  f,^(x)  and  f0(x)  to  the  sulfur  content;  and  f^(x)  and 
fj^Q(x)  to  the  volatility  index.  These  relationships  are  all  defined 
in  such,  a  way  that  the  particular  specification  will  be  satisfied  if, 
and  only  If,  the  function  f^(x)  >0.  To  sloplify  notation,  the  constants 
in  the  equations  do  not  refer  to  the  absolute  octane  number,  sulfur  con¬ 
tent,  or  volatility  index  of  a  particular  raw  gasoline,  but  rather  to 
the  differences  between  these  vsilues  and  those  of  the  specifications 
for  a  particular  motor  fuel  product. 

In  the  case  of  the  octane  ratings,  for  exsjqile,  equation  (V.l) 

Indicates  the  relationship  between  the  octane  nuak>er  of  a  particular 
blending  stock,  and  the  IXL  concentration  in  a  finished  product.  A 
new  function,  defined  in  order  to  indicate  the  nuia>er  of 

octane  points  by  which  stock  J  exceeds  the  speclfi eat ions  for  product 
k,  when  ml.  of  lead  are  added  to  stock  J.  'Rie  variable  t^^  is 
evidently  a  function  or  x^,  and  is  determined  throu^  the  following 


relation: 
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(V.8) 


‘J  ♦  ‘’J  *2J  - 


'J 


-  N. 


(k  -  1.2) 


l*x 


2k 


The  terms  Sj.  bj.  and  represent  the  lead  eusceptlbllity  constants 
for  stock  and  correspond  to  the  three  constants  in  eq\iation  (V.l). 

The  paraiaeter  stands  for  the  mlninuro  octane  number  requirement  con¬ 
nected  vith  gasoline  blend  k. 

In  order  to  ensure  meeting  octane  number  specifications  on  the 
premium  grade  gasoline,  the  following  condition  must  then  apply: 

(».9)  .5 


UM  -  J'-jiClj-Xj) 


)  J  -  41.  ?li 

)  60-66,  70 

)  i-  k5,  46, 
55s 


aimilarly,  to  make  sure  that  the  regular  grade  gasoline  blend  will 


pass  its  octane  number  specification,  we  must  have: 

(V.IO)  r^(x)  ♦  X3^  )  J  =  31,J6, 

-  ^,7*3^  -  0  ■  > 

)  80 

r 

% 

The  fSietlona  in  braobeta  that  follow  tlw  ooriahlos  <3^  ond  art 
Inserted  In  order  to  reflect  the  fact  that  with  these  two  gasollmaa, 

as  the  cutting  temperature  changes,  both  the  quantity  and  the  octane 

number  of  the  motor  fuel  component  will  be  affected.  (See  Figure  12. b.) 

This  feature  evidently  leads  to  negative  square  terms  involving  x^^. 

The  same  two  stocks  give  rise  to  quadratic  expressions  in  connection 

with  the  sulfur  and  volatility  specifications  for  regular  grade  gasoline. 

but  since  these  c  )aponents  were  assumed  to  be  absent  from  the  premium 

gasoline  blend,  those  equations  remain  completely  linear.  For  the  sulfur 

condition  reads: 

specification  on  the  premium  gvads  gasoline,  the 
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(V.u) 


fj(*) 


V‘°J  ‘'j’  ° 

j  -  41,  51,  60-66,  70 
k  -  45,  46,  55,  56 


'nie  constact  term  ^  represents  the  algebraic  excess  of  tiie  iwuclim 

allovable  sulfur  content  of  760O  gasoline  ever  the  anount  actusLlly  in 

associated 


stock J  .  "niere  is  a  similar  set  of  constants^ 


vlth  the  regular 


grade  gasoline  —  d  the  difference  between  the  allowable  sulfur  con- 
tent  In  the  76  grade  pxx)duct  and  the  content  In  stocky  .  n>e  foUcvlng 
relationship  controls  the  sulfUr  level  In  the  finished  regular  grade 
product) 

(V.12)  Cg(x)  .  E  -  X,^  ^  ♦  ^3k^J 

■  b<=  ° 

J  -  31,  36,  41-44,  51-54,  60-66,  80 


The  bracketed  terna  foUovlag  and  x^^  arise  out  of  the  cutting 
teaperature  relationship.  (See  Figure  12.  e.)  Rote  that  orar  the 
entire  range  between  and  the  terns  Involving  x^  bring  about  a 
strictly  concave  function.  'Riat  Is,  the  foUovlng  partial  second  de¬ 


rivative  is  negative  over  this  particular  Interval: 

'34*31 


(f.13)  ^  fe^x) 

T 


“33  ■  “■ 


*31 


Tbm  functions  f^(x)  and  do  the  Job,  respectively,  of 

governing  the  volatility  of  the  prenium  and  regular  grade  gasoline: 


f9(*) 


J 


'JL  ‘^3  -  -j'  *  f  ' 


kl^  - 


u 


J  -  hi,  H,  60-66,  70 


k  -  45,  4'  ,  36 


(V.r,)  fj_Q(x) 


■  ®32*31  ■  *37*36  -  ^ 


J  -  3I»  36,  sl-->,  bO-66,  80 


Here  a^in,  because  of  the  cutting  temperature -quality  relaticnehip , 
there  are  square  tema  connected  vith  and  x^^. 

The  gasoline  blending  model  has  now  been  reduced  to  a  constralned- 
aaxlffluin  ■athenatlcal  form.  The  problem  consists  of  choosing  a  vectc^r  (x) 
that  will  naalniae  g(x),  subject  to  the  aide  conditions  tliat  ’ 

and  also  subject  to  certain  upper  and  lover  limlta  on  the  Individual 
coerponeats ,  x  j.  The  payoff  function  g(x)  Is  given  by  equation  (V.U). 

The  upper  and  lower  lialtB  on  Individual  variables  are  stated  in  (V.b), 
and  the  10  restraining  inequalities,  f^!x),  In  (V.  ),  (V.^-  (V.U  ), 
(V.U),  (V.12),  (V.lU),  and  (V.is). 


In  this  fom,  the  problem  almost  fits  the  KUlm-Tucker  conditions. 

Both  the  payoff  function  and  the  restraints  are  concave  in  each  variable 

IH 

taken  separately.  furthemjore,  eight  the  ten  restraints  are  concave 
over  the  whole  s^ace  --  fj^(x)  .  .  .  f|j{x),  and  r,(x)  .  .  .  f^^(x). 

1  ft 

Kuhn  and  Tucker  give  the  foil. owing  definition  of  concavity:  "A  fmctlon 
la  coo  cave  If  linear  Interpolatioo  between  ita  values  at  any  two  points 
of  definition  yields  q  value  not  greater  than  its  actual  vslue  at  the 
point  of  Interpolation".  2^  clt .  ,  p.  48l,  The  function  is  said  to  be 
concave  In  each  variable  separately  if  the  two  points  of  definition  are 
always  chosen  so  that  only  one  f  the  coraponerts  of  the  vector  (x)  la 
altered.  A  function  may  be  concave  In  eacl;  variable  separately  v! thout 
being  concave  over  the  whole  space  of  points  (x). 


The  payoff  function  r(x)  and  the  two  octane  r.uir.ber  restraints  — 

and  f^(x)  —  do  not  strictly  satisfy  the  full  concavity  properties. 

It  WPS  decided,  nevertheless,  to  perform  the  computations  as  if  these 

also  satisfied  the  rinorous  Kuhn-Tucker  conditions.  Under  this  approach 

will 

there  is  always  the  danger  that  a  local  optimum  not  turn  out  to 

be  the  maximum  maxlmorum.  A  number  of  alternative  computing  schemes 

were  examine*!,  out  all  of  these  suffered  from  precisely  the  same  defect. 

The  primary  Justification  for  going  ahead  in  this  non-rigorous  fashion 

was  that  the  computational  results  could  always  t>e  checked  against  the 

experience  of  the  refinery  engineers.  It  was  felt  that  these  individual 
certainly 

he  able  to  detect  any  errors  of  a  gross  nature. 


would 
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U»  ”|*^Tic4l  analysis 

Outline  of  coMxitationfcl  proceduf 

Following  the  eesuaqilion  thet  this  gasoline  blending  model 

fits  into  the  foR^t  of  the  Kuhn-Tucker  theorem,  it  is  necessary'  to 

define  the  Lagrangean  function  0(x,u): 

10 

(V.16)  (^(x,u)  -  g(x)  4  r.  u^fi(x) 

i«l 

The  const ra ined-maxinum  problem  will  be  solred  if,  and  only  if, 
there  is  a  saddlepoint  solution  to  the  function,  0(x,u).  The  compu¬ 
tational  procedure  is  an  Iterative  one  -mt  each  step  t,  converting 
a  vector  x(t),  u(t)  into  a  new  vector  x(t  4  1}  and  u(t  4  1).  The 
solutions  observed  have  all  tended  toward  a  saddlepoint,  but  I  can 
give  no  strict  proof  of  the  necessity  of  this  convergence.  The  new 

vector  generated  is  never  exactly  "efficient**,  nor  is  it  *'attainable** 

19 

in  the  sense  of  Koopeans.  The  payoff  does  not  increase  monotonically 

as  in  each  successive  step  of  0antsig*s  simplex  procedure  for  linear 
20 

prograsning.  Despite  these  apparent  shortcomings,  the  algorithm  has 
given  useful  answers  in  the  cases  examined  to  date. 

The  solution  anist  be  started  off  from  some  initial  point  x(0},  u(0). 

In  principle,  this  m^  be  any  arbitrary  non-negative  vector.  In  practice, 
though,  it  is  possible  to  effect  a  considerable  reduction  in  computing 
time  if  a  good  initial  position  is  selected.  At  each  step,  first  the 
(u)  vector  is  detensined,  and  then  the  (x)  vector.  For  the  former,  the 

I9I  T.C.  Koopsians,  "Analysis ~of  Production  as  an  Efficient  Combination  of 
Activities",  Chapter  III,  Activity  Analysis  of  Production  and  Allocatiop.  p.  79. 

20.  G.B.  Dantsig,  *9(aximisation  of  a  Linear  Function  of  Variables  Subject 
to  Linear  Inequalities",  ibid^^ Chapter  XXI,  339'347* 
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basic  iteration  consists  of  tvo  steps,  end  for  the  latter,  of  four 
steps.  Tlie  problca  is  one  of  determining  ^  s  I)  -  t^t),  and 

A  Xj  f  Xj(t  1)  -  procedure  appears  cuaberscoe,  but  only 

six  minutes  of  CPC  aarhinr  time  are  actually  required  for  generating 
a  whole  set  of  A  u.  and  A  x  .  (The  required  total  number  of  steps  has 

«l 

varied  between  50  and  80.)  For  the  A  u^ ,  the  prooedxu«  Is  as  foUowat 
>(^|t|  >  0,  than  the  "oandldate”  A  •  -1^ 


1.  If  ^ 

a^i' 


If 


^Ct)  <  0, 

rtt} 


then  the  candidate  A 


note:  is  an  arbitrary  positive  constant. 


2.  If  uj't)  ^  A  <0,  than  A  *  “U^(t) 

If  u  (t)  +  A  u  >0,  then  A  u  =  A  u  (t) 

i  1  "  1  i 


For  the  A  x^  ,  the  procedure  beccoes: 
1 .  If ,  for  any  1 ,  fi  ^  t )  <  0  j 


21 


r  ^  «’•! 

<  0,  then 


=  0;  proceed  to  evaluate  A  x  j^.  Otherwise,  proceed 


to  step  2. 

2. 


If  0,  then  the  "candidate' 


^  "i 


^  "candidate"  A  *=  -k 


Mote:  kj  Is  an  arbitrary  positive  constant. 


21  .'^c  individual  A  xj  ai-c  determined  in  ascending  order  of  their  respective 
indices,  J_.  To  simplify  the  notation,  A  xj  .  is  used  to  indicate  the  next 
A  X.  that  is  to  be  determined.  Per  exsoplc^V  x..  follows  A  x_  in  sequence, 
and^for  J  »  31»  i®  Indicated  by  A  Xj  j^.  ^  ^ 


,-n 


If  J  *ad  If  Xj(t)  -♦•  A  Xj  ^  fJj. 


A  X 


J  J' 


-  ^j'.  proceed  to  «vmlu*te  A  x 


>1‘ 


If  J  ■=  31  or  36,  and  if  x^'t)  ♦  A  Xj  _•  .jj,  then  proceed 


to  step  U 


b.  If  ^  ^  31  or  36,  and  If  Xj’t)  ♦  A  Xj  •  0,  then 
Ax  *  -X  (t);  proceed  to  evaluate  A 

w  w  O 

.  \ 

If  J  »i  31  OJ"  3^),  fcnd  tf  Xj(t)  +  A  Xj  >0,  then  proceed 
to  step  U .  ^ 

a.  If  ^  «  42,  43,  43,  46,  ‘^2,  33,  53  or  proceed  to 

evaluate  A  x, 


‘J 


A 

Z  ^X^i 


-Av 


If  j  =  21  or  22,  calculate  /ToO  -  x.(t)  -  A  x  7. 

<3  ^  J 


If  this  eiqpreaalon  is  non -negative,  A  Xj  -Ax 


y 


Othervlse  A  x 


j 


3.00 


X  (t).  Proceed  next  to 
«/ 


evaluate  A  x 


•3  60 ,  3 . , 

ir  j  =  41,  44,  31,  34,  or><  .  .  .  66,  c«jlculate 


-  Xj(t)  -  A  this  expression  io  non- 

negative,  A  Xj  ^  A  Xj  .  Othervlse  A  Xj  ^  j  -  Xj't)  . 
Proceed  next  to  evaluate  Ax... 

for  the  A  .  the  explanation  of  this  ritual  is  stralghtforvard . 
Step  (1)  tells  us  to  decrease  u^  by  an  arbitrary  anount  if,  for  a 

"aisaLl"  change  in  alone ,  the  effect  will  be  to  decrease  0(x,u). 
Siallarly ,  there  Is  to  be  an  increase  in  u,  if,  for  a  gnell  increase, 
^  would  decrease.  Step  (2)  prevents  from  bccocilng  negative. 


rou;.  i.ibc'  .t. . 


F' r  tno  Xj  ,  Ui ;  .luctl  f  Icatl  m  In  a  rc/^ 


(Steps  (3) 

Rju';  (•♦)  are  the  obvifTOS  oncn  --  reopectivel^' ,  I'^wer  and  upper  Lljalts 
jn  the  InJ.vlduJLl  :<j  dictated  by  (V-ti).)  As  with  i.he  A  tne  sign 


is  prljaarLIi'  the  criterion  tliat  deteraines  whether  t''  tSLite 


a  rrontive  <  r  a  negati’/e  step.  Pi'clinlnary  snaiLl-acale  cadculAtljn* 

Irkllcated,  hcvever,  that  the  unq’ial  i  fled  "direction  of  ascent  (descent) 

rule"  would  lead  to  najnr  osctllatlona  ratlujr  than  t^  convergence  at 

equlllbrlun  vTitlaes.  For  that  reacon .  rule  (l)^  a  relajcat Ion -type 

principle, was  Inserted.  Tranelatod  Into  English,  the  stateaent  reada, 

"If  -ne  of  the  restraints  is  being  violated,  and  If  the  dlrectton-of - 

ascent  rule  would  lead  to  fin  even  greater  violation  of  this  restraint, 

then  the  particular  Xj  sliould  not  be  chsuge^^ .  *  The  relaxation  rule  Is 

.oarily  eiegsust  from  tlie  vlew^*'  int  of  pure  matheraatlca ,  azxi  there  are 

undoubtedly  other  possibilities  for  achieving  the  saac  effect. 

At  the  present  time,  leepite  their  plausibility,  there  It  a  serlmm 

dllerma  concerning  the  logic  of  ttiese  conqxitatlonal  techniques.  On  the 

on»>  hantl ,  I  have  been  'mable  to  prove  the  necessity  of  convergence  for 

either  tlie  rei-axat I 'r.  or  the  unaodifled  procedure.  On  the  other,  I  have 

CO  oroxam;  .9 

been  unable  to  constrict  m  for  a  case  where  the  k,^  and 

the  k,  any  be  made  arbitrarll:,'  laaall  .  Readily  granting  these  apparent 

j 

objections,  t.he  fact  rviminB  that  the  relaxation  teciuilque  has  givna 
useful  reaultB  in  e  eeriee  of  conpnitLng  .•"ons 
■•.1?  Dotailo  '  f  the  nurse  r  i  .aI  analy  ' 

Tl.ree  sets  of  cal cuLat  1  ons  have  been  perfonupd .  The  first  is 


to  bo  considered  the  basic  case  --  taking  as  paraaaeters  the  1953 

fhe  reader  should  note  tliat  ^  Is  being  nlnljalied  with  reapect  to  (u) , 
and  mxlinlzod  with  respect  to  (x). 
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refinery  rwtll lAt long,  product  sped fiCAt long,  gnd  rgw  aatgrial  avail- 

abllltleg.  The  socor^d  run  was  concerned  with  a  radical  change  In  the 

pdce  gtructuro  —  an  11.6/C  drop  In  the  value  of  prendum  grade  gasoline. 

In  the  third^  the  price  gtructure  wae  taken  to  be  the  aasM  as  In 

the  basic  case,  but  a  2.5  octane  point  increase  was  Imposed  upon  the 

conditions, 

preailujr  product.  For  each  of  these  departures  from  the  basics  the 

optimal  blending  schedule  becomes  considerably  altered. 

In  order  to  convey  some  Indication  of  the  computational  technique 
employed  here.  Figures  13. a  -  13. e.  chart  the  values  of  certain  of  the 
variables  during  the  last  twenty  steps  of  Run  #1,  i.e.,  between  t  40  and 
t  -  60.  According  to  Figure  13. a,  the  payoff  and  the  aggregate  make  of 
premium  and  regular  grade  gasoline  had  become  etabillsed,  and  "rippled" 
stsadlly  sbout  certain  equilibrium  levels.  Virtually  all  of  the  inde¬ 
pendent  variables  exhibit  similar  rippling  around  an  squill brl’oin  point, 
but  there  were  exceptions.  The  worst  drifting  occurred  in  the  case  of 
respectively,  ths  quantities  of  stock  #42  in  regular  and 
premium  grade  gasoline.  Their  time  paths  are  shown  In  Fiinire 
Ths  more  typical  pattern  was  exhibited  by  the  two  TEL  concentration 
variables,  and  x  ,  and  is  reproduced  In  Figure  13. c.  Two  of  the 
Lagrangean  multipliers,  u,  and  u  ,  connected  with  the  two  octane  number 
specifications,  are  shown  In  Figure  13. d.  The  octane  number  functions 
themselves,  fc(x)  and  f,(x),  are  given  In  Figure  13. e. 

Before  Cutting  off  any  Individual  computation  run,  it  was  possible 
to  perform  an  additional  check  upon  the  stability  of  the  solution. 
According  to  the  Kuhn-Tucker  theorem.  In  onler  for  tke  vector  (x)°  to 
be  a  saddle  value  solution,  the  following  condition  must  hold  for  all 


23.  Og.  cit. ,  p.  482 


''HC 
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(V.l,’) 


-L  0.  «j  if. 


0. 


In  the  present  cs»e,  where  the  Xj  are  bcmnded  froa  above  as  wall 

24 


ftc  below,  this  condition  on  the  equilibrium  becooes:  ^ 

(V.l3)  a.  If 

b.  If 

"sing  ronditlcn  (V.l8),  it  Is  pocolble  to  yerfora  a  rough 


<  0,  then  X  “  .  0,  and  x°  > 

J  5  T  •=  J  ' 

'*  u,  then  x^  ■  (^ 


check  on  the  pr>-clsloc  cf  the  solution.  For  those  x  that  are  at  their 

J 

lower  I'nltB,  o  ^  _  G.  For  these  that  are  at  their  upper  limits, 


^  X 


f 


j 


>  0 .  And  for  those  that  are  neither  at  their  upper  nor  at  '.hei' 


’  X 


Icwer  Halts,  ^  *  0.  In  this  latter  case,  it  is  not  to  be  expected 


IX 


J 


that  *hc  partial  derivative  vUl  be  literally  at  a  rero  level.  Rather, 
since  this  is  an  approximative  metn  >d  of  Bolution  i-vi»  derl'.'s.tlvc  should 
be  In  the  oeigliborhood  of  tero. 

Table  V.2  cc^talns  the  actual  values  for  the  xj(6o)  and  for  the 
normalised  values  of  the  partial  derivatives  6o) 

variables  that  art  at  their  lower  limits 
as  Is  to  be  exjectei. 

negative  partial  derivatives, A  On  the  other  hand,  the  single  variable 


All  of  those 

exhibit 


e  case  of  (V.l3a),  where  ^  -  SI  or  36,  the  lower  llailti  q?  Kist 
erpreted  as  equivalent  to  aero.  For  {V.l8b),  where  «  21  or  22, 


In  the 
be  Interpi 

the  upper  limits  of  3. 00  ml.  are  to  be  regarded  as  the  appropriate  (x  . 

For  those  xj  on  which  there  Is  no  specific  Imposed  upper  limit  --  tMt  Is, 
42.  I3.  4‘^ .  46,  “^2,  63.  '^6,  Sind  ‘^6  --  the  conditio©  (Y.lSb)  does  not 
applG' . 
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Table  V.2  tomparlicm  of  Partial  Derivative*  for  t  6o,  Rrai  #1 


1 

Lcfwer 
Limit 
on  xj 

Xj(60) 

l»pper 

Limit 

on 

3^(60) 

:95tj(66) 

A<ijusted 

Optimal 

21 

0 

1 .  )4 

3-00 

ni'. 1 '  lor  in 

1.5711  of  co*t  of  1a 

prmniun.  gasoline 
1.92 

22 

0 

2.)e 

3 .00 

-0.6l%  of  co*t  of  ly 

2.96 

31 

l,90C 

1.900 

^050 

n.!. .  i  .iter/fa.  .o:.  ir. 

-22.54lt  of  realization  on  1  barrel  of 

regular  gaso.ine 

1,900 

36 

1,320 

1,320 

2,060 

prralUD  gafiollne 

-0 . 34  ■ 

1,320 

Ul 

0 

0 

i.ono 

-5 -5911  of  realliatlon  on  1  barrel  of 

0 

42 

0 

540 

2,910 

premium  gaeollne 

^.80 

620 

^♦3 

0 

400 

3,U6o 

-0.40  "  ■ 

380 

44 

0 

c 

1.790 

.0.62 

0 

45 

0 

2.3/) 

2,910 

2.35 

2, 2^X3 

46 

0 

350 
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1.09  * 

845 

51 

0 

0 

no 

-4.4llt  of  reallratloQ  on  1  bforel  of 

0 

52 

0 

0 

1,240 

premium  gasollx^ 

-3.41 

0 

53 

0 

200 

860 

-0.05  ■  ■ 

200 

54 

0 

4  0 

4'X) 

-0.06 
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55 

0 

100 

1,2U0 

0.26 
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56 

0 

oOO 

360 

1.21 

600 

6o 

0 

0 

6.150 

-b.901l  of  realltatloo  on  1  barrel  of 

0 

ol 

0 

0 

3,‘*30 

premium  gaaoline 

-5.90 

0 

62 

0 

0 

l,46o 

-9.02  • 

0 

6'3 

0 

0 

3,200 

-6.55 

0 

64 

0 

0 

1,600 

-^.10  ■  • 

0 

65 

0 

0 

11,230 

.<.'94 

0 

'"6 

0 

0 

7,440 

-7 . 35  ■  ■ 

0 

tliflt  at  Its  uppxr  limit,  x,  ^ .  does  n^t  have  a  partial  derivative 

of  the  proper  sign.  'Hilc  derivative  is  slightly  negative,  and  in  fact, 
during  the  sequence  leading  up  to  t  ^  60 ,  x,  ,^(t)  had  oscillated  betwaen 
UYO  B  ^CD  and  ^pK).  There  aopears  to  be  a  coincidence  between  the  upper 
limit  on  this  particular  variable,  ^nd  the  point  at  which  0  vaniahss. 

In  the  caae  of  those  variables  that  lie  neither  at  tiialr  upper 
nor  their  lower  limits,  the  naximunt  absolute  deviation  fron  lero  aocnints 
to  -f  the  value  of  one  barrel  ii'  j-renrl’in  grade  gasciine.  This 

Icv.nti^n  occurred  In  the  case  of  x  ,  the  variable  shewn  in  figure  V.lb 
ft  -  .  nad  been  slowly  drifting  upward.  It  did  not  seem  worthwhile  to 
continue  t.'.e  iterations  in  order  to  refine  ujon  this  one  estimate. 

T^e  blending  schedule  selected  at  itep  ou  is  neither  literally 

'’attair.^b  ,  e*' .  Iio 

■  .ent"  nor  ^  payoff  g(x)  Is  balnv  the 

a.-.:  for 

le  ^el  recorded  for  step  ^5*/  thla  reason,  the  final  blending 

schedule  was  taJeen  over  froo  that  earlier  step.  Coitain  minor  cor¬ 
rect  loos  were  made  by  hand  in  order  to  bring  about  a  feasible  and 
efficient  solution.  The  resulting  vector  is  Listed  in  the  final 
column  under  the  heading  "Adjusted  optimal  Xj". 


.  Results  of  the  numerlca.!  4nalj8i8 


5 . 1  Interpret4t ion  cf  results  —  bleniiry:  pro^.raJpa 

The  operations  sche<iules  Jerivel  frorr.  the  three  computation  runs 
7  reser.tel 

are/<  in  Taole  V.3.  For  reference  purposes,  the  first  three 

columns  Indicate  the  program,  that  was  effective  at  the  beginning  of  the 
year  1953.  The  other  sets  give  the  results  for  run  #1,  the  basic  case; 
run  #2,  a  iecrease  of  11.5'^i  in  the  realization  on  premium  grade  gasoline; 
ani  run  an  increase  of  2.*^  octane  numbers  in  the  specification  for  the 
premiuir.  proluct. 

The  outstanding  features  of  these  lifferent  sche  iules  are  suinraarize: 
in  bar  charts,  Figures  Li*. a.  -  L'*.!.  Of  all  four  situations,  thf*  outp^ut  of 
7b00  gasoline  is  Largest  in  the  case  of  run  #1,  nnd  smallest 

under  '’1953  operations".  Just  the  reverse  is  true  for  the  regular  graiie 
f^sol i  ne.  Its  production  is  largest  unoer  t.ie  19^3  operatlor-'  .lin,  ari 
smallest  in  the  basic  case.  Going  along  wit^  tnis  considerable  shift  Ir. 
the  general  pro<luction  strategy,  there  is  an  apparent Lr  small  shift  in  pay¬ 
off,  for  the  naslc  case  would  yield  a  return,  to  the  company  of  only 

25 

percentage  points  in  excess  of  that  shown  ir  actual  practice.  At 
first  glance  this  figure  of  1.99  loes  not  seen  impresalve,  but  neverthe¬ 
less,  luring  the  course  of  a  single  year’s  operation  a  differential  of  this 
magnitude  amounts  to  several  millions  of  dollarsi  From  the  comp^any’s 
standpoint,  it  would  look  eminently  worthwiille  to  shift  over  to  the 
blen>dinr  schedule  derive  1  from  run  #1. 

25.  In  order  to  avoid  revealing  proprietary  information,  no  absolute  dollar 
figures  sre  presented  here.  Both  the  unit  and  the  total  refinery 
realltatlons  are  stated  in  terms  of  index  numbers. 
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Unhappily,  the  life  3f  the  operatlona  aiMdyst  doei  not  run  a 
slaqiie  courae .  frctt  the  very  beginning  of  the  rtudy,  It  bad  been 
acknovLedged  that  the  conpany  could  not  realiatically  expect  to  aall 
uni 'ml  ted  quantities  of  prealua  gasoline  at  the  prevaU.lng  refinery 
realization.  At  the  same  tlae,  it  had  net  been  fully  appreciated 
how  large  an  increase  in  the  production  of  this  Itea  might  be  in 
order.  (T^  cutput  In  toe  one  case  is  almost  double  that  of  the 
other.)  For  this  reason,  it  seems  desirable  to  Incorporate  an 
additional  fomsl  restriction  into  the  problem,  and  to  require  that 
the  total  production  of  premium  product  remain  belcw  scam  upper  limit, 

WO'.  1  Oft 

.  In  this  way,  It  ,ponalble  to  detemine  how  large  an  increase  in 
is 

profits/  attainable /even  under  the  limitation  of  current  sales 

of  the  one  product.  B:-'  varying  this  ui^per  limit  parameter, 

i'  1  p:i’  'O'^.J 

trace  out  the  effects  of  altcmatlre  sales  volumes 

•-O’  I 

upnin  the  ^  payoff. 

Rather  t:ui.n  i-evlae  the  ooraputlng  setup  in  order  to  impose  a  direct 
marlcetlng  restriction  if  tnls  natu’*'',  It  ceened  r  .TH.re  eror.coil^'al 
alternative  to  calnlnte  a  second  point  along  the  empany's  "supply 
curve"  for  preni'jci  gasolia".  Tha*  Is.  for  run  f2  all  parameters  wore 
maintainod  at  *he  same  Level  as  in  run  #1.  exceot  for  the  unit  reall- 
zatl  )n  upon  ’'OC  gasoline.  Tills  was  dropped  to  66.4?^  of  its  previous 
level,  "he  -hange  had  tne  effect  of  reducing  th#*  dlffe.'^ntial  between 
prcial  ....  tun  I  .  »  g  lao  grade  gasoline  to  Less  thsus  half  of  Its  prevlcrus 
funount .  ?b‘.B  alteration  in  the  rrlce  relations  ha.-"  a  result  ansklogo^is 

t  -1  that  if  a  ulnect  o^iantlty  limitation  io<jn  xAvt  najcc  of  premlun  gaac - 
lint*,  ihe  out  Imal  output  oi’  this  nrnluct  was  re<.iuce<i  from  P  d* 
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down  to  26,36^,  while  the  total  of  regular  grade  76  gaaollne  went  up 
froa  13,670  B/CD  to  33,390.  In  this  second  case,  not  only  prenlun 
^isollne  hut  also  cutter  stock  production  declines.  It  now  becooes 
profitable  to  Incorporate  nost  of  the  hsawy  gasolines  Into  notor  fuel 
rather  than  to  blend  then  into  residual  fUsl  oil. 

The  results  of  run  #2  are  equivalent  to  those  that  would  be 

derived  by  esgtloying  an  outright  sales  limitation  of  26,3^3  B/CD.  'Rie 
Twrlce  tc  76CC  —  —  Luuicaies  me  aarglnai  oppor¬ 

tunity  cost  of  this  product.  In  other  words,  if  the  narketing  restriction 
were  raised  by  one  barrel  frost  26,38^  B/CD,  and  if  aaterials  were  avail¬ 
able  for  mking  this  extra  barrel,  the  company  would  be  able  to  obtain 
an  additional  payoff  of  8d.U2)(  of  the  current  average  realisation  on 
7600  gimoliae.  In  the  Jargon  of  econooists  this  value  is  termed  the 
"shadow  parioe*  of  the  premluB  product.  Evidently,  there  is  a  marked 
gsp  here  between  the  shadow  price  and  the  average  realimtion  upon  the 

26  -  -  -  .  •  - 

followers.  This  groop  hM  stresssd  ths  fast  BiMi  sudi  digersiaasl— 

sre  symptaaatic  of  aarket  restrainta.  At  the  same  time,  ^  thsm- 
selvee  these  discrepancies  constitute  nc  proof  of  either  the  absence 
or  the  presence  of  workable  coagietition.^  Prom  the  viewpoint  of 
internal  asnageaent,  they  do  provide  the  cai^mny  with  an  indication 

'^8ee  especially  aT  P.  Lemer,  ^Sie  Coneept  of  Monopoly  and  the  Msasure- 
nent  of  Mr^nopoly  Power",  Review  of  Economic  Studiea,  June  193*^,  I, 

pp.  157-175. 

7  ^ 

For  a  discussion  of  this  concept,  see  J.  M.  Qark,  toward  a  Oonoept 

of  Workable  Cos^tition",  Aneriean  Economic  Review,  June,  19^,  Vol.  XXZ, 
Ho.  2,  pp.  2UI-256. 
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of  the  groes  gains  that  night  be  aehieTed  bj  reorienting  sales  polio/. 

One  of  the  fsTorite  sTenues  for  sa3.es  rivalr/  in  the  sascline 

industry  consists  of  alteration  in  the  specifications  for  the  house- 

brand  product.  These  aodifi cations  occasional!/  take  the  fom  of 

changes  in  the  eolatilit/  index,  but  more  t/pieall/  center  about  the 

octane  ratings.  For  this  reason,  it  «as  of  particular  interest  to 

explore  the  effects  of  an  octane  number  increase  upon  the  optiaal 

the 

blending  schedule  and  upon^pa/off.  la  run  #3,  it  uas  assuaed 

v'nat.  me  specification  of  the  premium  gasoline  had  been  increased 
b/  ?.5  octane  points  —  a  plausible  Jump  in  this  requireasnt.  Keferring 
back  to  the  bar  charts.  Figures  14. a.  -  14. d.,  the  reader  can  compare  the 
results  with  the  basic  case:  a  decrease  in  total  pa/off,  a 

noderate  decrease  in  the  output  of  presdum  7600  fuel,and  an  increase 
in  the  output  of  76  gasoline.  The  make  of  cutter  stock  reaains  the 
same  magnitude  in  both  eases. 

The  oxitstanding  difference  between  this  situation  and  the  others 
consists  of  the  high  eth/1  fluid  consumption.  It  aitounts  to  6,755 
liters/da/  here  —  about  a  one-third  increase  as  against  5,026  liters/ 
da/  in  the  basic  ease.  The  IVL  concentration  leTol  is  direetl/  up  at 
the  prescribed  limit  of  3*00  ml. /gallon  for  the  premium  grade  product, 
but  despite  this  fact,  it  is  onl/  at  2.60  ml.  for  regular  gasoline. 

Like  the  results  of  the  reformer  anal/sis  in  chapter  IV,  these  compu¬ 
tations  do  not  /ield  much  comfort  to  those  who  would  put  faith  in  an 
oftem-quoted  rule  of  thumb  —  that  the  optimum  blending  polio/  is  to 
assign  the  maximua  allowable  TIL  concentration  to  each  gallon  of  product. 
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Only  in  rnnkinr;  prerium  ‘'raiie  f^solins  in  run  #3,  and  in  malting  the 
regular  product  in  run  #1  was  the  3*00  nd.  limitation  “pproachcd. 

In  all  other  cases,  it  paid  to  employ  less  than  this  concentration. 

Not  only  is  the  ethyl  fluid  consumption  closely  connecte<l  with 
the  octane  specification  of  the  two  motor  fuels.  It  is  also  sensi¬ 
tive  to  a  variation  In  the  price  structure  for  premium  grade  gasoline. 
Both  in  the  1953  operations  and  in  run  #2,  the  7600  product  commands 
a  low  shadow  price,  and  by  cotr.parison  with  the  basic  case, there  is 
nCvui  uxil;  ,1/  a  low  consumption  oi  etnyi  fiuld.  Kuns  fl  and  #3 

indicate  that  even  during  the  ordinary  peacetime  course  of 
market  conditions,  there  is  a  high  degree  of  flexibility  poaslble  in 
the  use  of  this  item.  During  wartime  emergencies,  it  follows  that 
there  are  substantial  opportunities  for  cutting  down  on  the  use  of 
the  material  in  motor  fuel  in  order  to  release  TEL  for  aviation  blends. 

An  interesting  sidelight  to  the  motorist:  Under  the  schedule 
set  up  in  the  basic  case,  there  is  actually  a  lower  TEL  concentration 
in  the  7600  "ethyl”  gasoline  than  in  the  regular  grade  product.  This 
does  not  mean  that  the  consumer  fails  to  get  his  money’s  wort-h  on  the 
premium,  product.  The  ethyl  fluid  concentration  by  Itself  is  meaning¬ 
less  as  an  indicator  of  gasoline  quality.  In  orier  to  determine  the 
knock  ratine,  it  is  also  necessary  to  specify  the  composition  of  the 
product  —  the  rexative  amount  of  components  lerived  from  straight-run 
gasoline,  from  catalytic  cracking,  from  tnermal  reforming,  and  so  on. 

5.2.  Interpretation  of  results  —  Larranr.ean  multipliers 
Because  of  the  need  for  withholdln"’  cost  and  revenue  figures, 
it  is  not  possible  to  Indicate  the  precise  nur..orlcal  values  for  the 


ten  Lagrangear.  muitij-liers,  u^.  ‘.evarthe-asa ,  the  ger.erai  nature  of  severa'. 
results  aay  M  meiitiotied.  Six  of  the  iaif;rangean  raultipiiers  tumeu  out  to  be 


poaltlw  iB  all  thrM  niaat  u^,  tha  aaWrlal  baLMce  for  atoek  flt$} 

Uj^,  the  aaterlal  balance  for  atock  #*j3;  u^,  the  octane  nunber  of  ^b00 

gaaollne;  u^,  the  octane  number  of  ^6  gaanllne;  u^,  the  volatility 

Index  of  76OO;  and  Uj^^,  the  volatility  index  of  Jf>.  Since  the  Ing- 

rac«ean  aultlpliera  were  alvaya  positive  In  these  six  casee,  the 

restraint  conditions  were  always  satisfied  exactly.  It  never  paid 

to  give  away  octane  nuabers  or  volatility  Index  points,  and  it  never 

paid  to  eaploy  atock  #k2  or  #53  aa  cutter.**" 

A  rather  different  general  result  held  true  for  u.^  and  U0  •-  the 

"ahadow  prices"  aeaoclated  with  the  sulfur  content  of  the  two  finished 

for 

gaeolines.  These  two  \i^  consistently  turned  out  to  be  zero^ 


.  vtah  tha  Uandiag  atocka  currsattly  swwiilabLa,  the  ecmfmoj  aa  a  abola 
aalatataad  a  eoarftarta^la  poaltlsa  wltk  raapact  to  aaatljag  ita  auLflar 
aped  float  ions.  Ttiis  reault  had  not  been  foreseen  In  advance,  and 
it  turned  out  to  be  fortunate  that  the  requirement  was  pbaaaed  In 
tarae  of  a  imxiiair  sulfUr  content  rather  than  an  exact  one.  Rad  It 
been  specified  as  an  enact  requirement,  a  dlatinctly  non-optlaal 
prograa  might  have  been  generated. 

Cooceming  the  two  material  balance  equations  associated  with 
stocks  #h3  and  #52,  no  general  rule  can  be  laid  down.  Both  u^  and 
turned  out  to  have  sero  values  in  runs  #1  and  #3,  but  achlavad 
20 - 

Tkble  V.3  shows  certain  quantities  of 
stock  #53'^^*^  as  nitter  in  the  case  of  rune  #2  and  #3.  This  is  noo' 
optimal,  buC* is  s  natural  consequence  of  the  approximative  type  of 
aolutloe. 
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pcGitive  levels  In  the  case  of  ran  Tliat  is,  for  sm  optinial 

SvJ.ution  under  the  parasKters  of  the  second  run,  uelUiex-  ux  these 
stocks  should  be  anployed  as  cutter  naterlal. 

The  uj  are  of  interest,  not  only  in  indicating  whether  or  not 
the  restraint  conditions  were  satisfied  ao  exact  eqitallties,  but 
also  ae  a  nethod  for  evaluating  the  use  of  ad^lltional  blending 
oaterial  in  either  of  the  two  motor  fuel  end  products.  These  vari> 
ahlcs  establish  a  scale  for  deteraining  the  dollar  bonus  or  the 
penalty  that  should  be  attached  to  small  incremental  quantities  of 
a  material  that  either  exceedf  or  i'aJ  .s  chort  of  the  finished  pro¬ 
duct  specifications.  The  first  use  of  these  Lagrangean  nultipliere 
was  to  test  the  earlier  assuaptions  as  to  which  blending  stocks  could 
be  excluded  froo  consideration  as  cociponents  for  7aOC  and  which  others 
could  be  rxxled  out  as  coraponents  in  the  7C  gasoline  product.  Under 
the  existing  configuration  of  blending  stocks  and  market  coetditions, 
these  assinaptluns  did  <±'Ck  out. 

It  is  e]q>ected  that  the  Lagrangeans  will  also  be  of  use  daring 

the  interval  before  the  coopany  is  in  a  position  to  schedule  the 

blending  problem  simultaneously  with  the  intra-ref luery  cracking  and 

converclon  operations.  The  ccopany's  Los  Angeles  refinery,  for 

cxanple,  has  a  polyfonaer  unit  that  has  econosiic  features  siailar 

to  the  thenaal  reformer  discussed  in  Chapter  IV.  As  with  the  Kellogg 

Coi^iany's  refonaer,  there  is  a  continuing  problem  of  selecting  tbs 

optinua  operating  coodltSons.  The  higher  the  reactor  temperature, 

the  better  will  be  the  octane  reting  of  the  polyforaer  gasoline, 

-ver  .juAJ’.tlt-  ,  w  . 

but  the^  will  be  they)  The  outstanding  difference  between 
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th«  two  problems  is  the  fsct  that  the  one  reformer  was  coneidered  as 

an  isolated  piece  of  conversion  equipment ,  and  that  the  other  forms 

just  a  single  component  within  an  integrated  company.  In  this  second 

case,  because  of  the  possibilities  for  blending,  there  is  no  simple 

linkage  between  the  market  realisation  on  the  end  products  and  the 

operations  of  the  one  unit.  Ideally,  the  company-wide  optisdsation 

doubtful  whether  this 

should  be  treated  as  a  single  calculation,  but  it  is  a 

will  be  achieved 

/fwithin  the  foreseeable  fdture.  Palling 


wf  t'lila  UmI,  lb  •wulvi  be  a  uBeiul  svep  lo  aLvaca 

the  polyformer  unit  problem  through  the  shadow  prices  derived  from  the 

blending  ealcxilation.  Theory  assures  us  that  this  is  a  legitimate 

assumption  if  the  company,  as  a  whole,  is  at  an  optimal  point  in  a 

concave  progranmir^  problem.  Theory  does  not  assure  us  of  this  result 

if  the  company  fails  to  be  in  equilibrium,  but  it  strongly  suggests 

that  the  same  principle  holds  true  even  when  the  oonpany*s  current 

29 


position  is  not  an  equilibrium  one. 

5.3  Company  policy  implications 

This  much  is  already  evident  from  the  computations:  that  it 
have  paid 

would  A^or  Union  Oil  to  forego  several  siilllons  per  annum  if  a  sales 
have  bean  for  disposing 

strategy  could  /Revised  ^of  an 

increased  volume  of  premium  gasoline.  For  this  reason,  the  problem  was 


29.  For  a  discuseion  of  the  potentialities  for  decentralised  decision¬ 
making  through  ths  use  of  shadow  prices,  see  Oskar  Lange,  On  the 
Iconomic  Theory  of  Socialism.  Minneapolis,  193S«  end  also  A,  P.  Lamsr, 
Bconoadcs  of  Control.  Now  fork,  19U».  Also  see  C.  J.  Hitch,  "Sa^ 
optimisation  in  Operations  fhroblsms".  Journal  of  tha  Ops rations 
Resaarch  Socisty  of  America.  Vol.  I,  No.  3,  May,  1953,  pp.  P'’-99| 
and  A.  W.  Marahall,  "A  Mathamatlcal  Not#  on  Sub-optlmlsatlon",  Ibid. . 

pp.  100-102. 
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takan  up  with  the  company**  sale*  senrice  tnanatTer.  A  number  of 

alternatives  were  discussed  with  him:  price  reductions,  quality 

Improvements,  and  an  advertising’  campaign.  He  agreed  that  all  of 

these  steps  could  conceivably  increase  the  sales  of  premium  grade 

gasoline,  and  that  the  potential  increase  in  revenues  might  exceed 

the  extra  costs  that  were  entailed.  He  was  pjartlcularly  concerned 

of  forecasting 

with  the  Ilf ficultie%^ 

the  response  of  consxuoers  to  such  changes  in  strategy, 
but  did  not  bring  up  the  parallel  question  of  how  the  company*s 
innedlate  rivals  might  react. 

Rather  than  risk  a  major  departure  in  price  policy,  he  indicated 

Union  Oil 

that  it  would  be  useful  for/]  to  explore  the  possibility 

of  altering  its  retailer  gasoline  discount  structure.  On  any  sales 
in  excess  of  certain  goals,  the  service  station  operators  might  be 
allowed  an  additional  discount  from  the  tank  wagon  price.  By  read¬ 
justing  the  pro^iuct  quotas  and  the  sales  bonuses,  he  believed  that 
it  would  be  feasible  for  the  retail  operators  to  shift  sales  in  the 
direction  of  7600  and  away  from  76  gasoline.  This  policy  is  one 
that  need  not  iimnediatsiy  provoke  retaliation  from  the  company's  major 
-iva.s,  and  it  does  not  depend  upon  a  large-scale  publicity  campaign. 
There  is  the  added  advantage  that  the  strategy  can  be  altered  at 
short  notice  without  making  consumers  aware  of  a  delltarat*  attempt 
to  influence  their  choices. 

Even  aside  from  major  departures  in  company-wide  sales  policy, 

was 

the  gasoline  blending  ncdel  yielded  on*  result  that/^  expected  to 
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produce  ea  lADadlats  reduction  in  cost.  In  nil  throe  conputing  runs, 

It  turned  out  to  be  optlsnl  to  cut  motor  guooline  coiqponent  #31 

i.e. . 

^  C  O 

the  lowest  possible  end  pointy  at  300  mthcr  then  ilie  current  332  F. 
This  chsnge  in  operating  xuactioes  could  be  effected  vithout  mklng 
any  radical  alteration  in  the  end-itee  product  mix.  At  the  aeoe 
tine,  the  isQ^roveaent  in  payoff  that  would  result  from  this  single 
step  amounts  to  a  suq  which,  in  less  than  two  days,  could  Justifr- 
the  entire  computing  machine  expense  involved  in  the  study. 
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6.  gumaniry  mui  conclualons 

grT^  Um  Msuapti>>a  «vf4fixet'  availabilities  of  rav  blendijag 

1 

ziatcriaic  aal  of  fixed  realizations  upon  flnislied  products,  the  Union 
Oil  Conqpony's  gasoline  blending  operations  were  formulated  in  tenns 
of  B  constrained -aaxiniur.  problem. 

Tfie  outctaading  result  of  the  noiaerical  eauLlysls  was  ^he  #Be» 

Uiiat  it  would  be  worth  several  millions  of  dollsirs  i>er  annua  if  the 
company  were  able  to  shift  sales  so  as  to  sell  additional  quantities 
of  prealua  grade  gasoline,  a  corollary  of  this,  It  foUovs  that 
the  cooiany  ban  consldcirable  latltatie  in  Its  ability  to  shift  production 
between  the  two  grades  of  la^tcr  fuel.  MarlQeting  problems  aside,  the 
Ideal  pattern  of  production  is  quite  sensitive  to  the  price  structure 
that  provaLLs. 

A  st'cond  result  of  the  calculations  •-■as  the  finding  that  the 
cmapan:,'  ccfuld  increase  the  octane  specification  of  Its  premium  pro¬ 
duct,  and  that  even  aft.er  such  a  change,  the  optlmsl  production  rate 
of  this  one  item  would  rsEmla  in  excess  of  cunwit  sales. 

Rsithsr  of  these  results  leads  to  tbs  laplloatlom  that  the  on^any 
■iTugbt"  to  change  Its  price  policy  or  its  product  specifications. 

Ml, lor  changes  of  this  natxire  cannot  be  evaluated  unless  the  coaomny's 
executives  also  take  some  account  of  the  response  of  consuasrs  and  of 
other  petroleum  refiners.  Readily  granting  the  uncertainty  that  must 
attend  any  such  conjectures,  the  fonoal  oathematical  treatment  seems  like 
a  proailslng  approach  to  the  blending  problem.  It  did  not  turn  up  vith 
solutions  that  are  absolutely  unobtainable  by  the  custonary  trial -and- 
error  techniques.  Rather,  it  brougiit  to  li^t  nuaermis  assuiqptionB 
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that  'uul  previously  Iain  hidden,  and  holds  out  the  possibility  for 
cl:>sei  x>nilxi&tlon  between  the  manufaicturlng  operations  sod  the 
activities  of  other  bxv&ches  of  the  corporation.  It  Is  particularly 
Isportent  for  the  aales  group  to  have  c]uantltatlve  inforaatlon  avail¬ 
able  on  the  possibilities  for  shifting  between  various  classes  of 
refined  products,  and  for  changing  the  specifications  of  those  products. 
It  is  equall>'  significant  for  them  to  know  the  cost  implications  of 
such  clAnges. 

Sxperiencs  within  the  iJhlan  Oil  Cosipsu^  suggests  that  aatheaatloal 
SK^els  have  considerable  utility  in  dealing  with  conplex  production 
probletas.  ’Rils  study  also  indicates  that  ecoocolsts  and  businessmen 
alllfe  are  a  long  distance  from  having  any  satisfactory  theory  of  how 
a  mrket  vlU  respond  to  s  change  in  the  oviBpmay't  sales  polieles. 

In  the  absence  of  such  teovledge.  It  Is  dangerous  to  place  too  nuch 
reliance  on  f or  h1  technlq[ues.  This  saase  danger  exists,  though,  with 
back-cf- the -envelope  calculations.  In  none  of  these  natters  can  a 
aathenatleal  nodel  make  important  business  decisions  by  Itself.  Iha 
Boet  that  can  be  claimed  is  that  it  will  suia»riae  the  volune  of 
Inforaatlon  that  executives  have  to  act  upon. 


